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1.  Preliminaries. 

Let  V'  be  a  real  Hilbert  space  and  n  :  V  x  V  — >  R  a  continuous  bilinear  form.  It  is  a 
more  or  less  standard  result  that  if  n  satisfies  the  conditions  .  Ji 


>"X  s,,P’ir~iT  7T37  >  c  >  0.  |  J 
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where  c  is  a  constant,  and 

(1.16)  (n(f.tr)  =  0.  Vr  6  1'}  te  =  0 

then  there  is  L  €  GL(V).  unique,  such  that 
(1.2)  a(r,  u>)  =  (Lv.  u>), 


By . 

Di.t  ib  i,„-  / 


Av.uidbiMy  Cones 


Dist 


m. 


Aviti!  j  'djor 

Special 


and  hence  given  (  €  V",  there  is  a  unique  solution  u  6  V  of  the  problem 


(1.3) 


«(?/.  r)  s=  f(r),  V  r  €  \ 


Indeed,  existence  of  L  C  £(V)  in  (1.2)  follows  from  the  Riesz  representation  theorem, 
and  (11a)  easily  implies  that  L  is  one-to-one  with  closed  range.  Next,  from  rlosedness  of 
rge  L  and  the  relation  rge  L  =  (ker  L *  )L.  and  since  condition  (1.1b)  implies  kerl*  =  {0}, 
we  infer  that  L  is  onto  V'.  i.e.  L  6  GL(V)  by  Banach's  theorem.  Next,  if  y  represents  (  via 
the  Riesz  representation  theorem,  that  is.  f(r)  =  {g.  r),  Vr  €  V.  »/  =  Z.-,y  is  the  unique 
solution  of  (1.3). 


I 


Conditions  (l.la/b)  include,  but  are  not  limited  to,  the  rase  when  a  satisfies  the  hy¬ 
pothesis  of  the  Lax-Milgram  lemma 

From  now  on,  we  assume  that  V'  =  where  Q  C  Rn  is  a  bounded  open  subset 

with  Lipschitz  continuous  boundary  I\ 

Let  {Qh  )i<K<N  be  a  collection  of  open  subsets  of  Q  with  Lipschitz  continuous  bound- 

.v  _ 

aries  TV,  and  such  that  O  =  0  if  A'  /  L  and  V  ft*  =  ft  For  convenience,  we 

K—  1 

shall  refer  to  ((Ik)  as  a  "partition"  V  of  fl. 

We  shall  denote  by  J  the  set  of  pairs  ( K.  L )  with  I\  /  L  such  that  meas  (T/,  0  Ti )  >0. 
Here,  “meas"  refers  to  the  r p  .  or  equivalent ly  r p-  Lcbesque  measure.  Also,  we  set 


Ik  =  {1  <  L  <  .V  :  (AM)  e  £},  1  <  A*  <  ,V, 


and 


rat.  =  rK  nrL.  v(K,i)€i. 


Let  Tr  he  the  disjoint  union 

r*=  U  r,w. 

(AM.  167 

(hence,  in  TjShL  and  T lk  are  different  although  they  coincide  as  subsets  of  Rn). 

We  introduce  the  space 

t-'ti'i)  =  n  l'( rM.». 

lit  141 

Elements  of  L2(Ti)  thus  are  collections  A  =  (A kl)  with  A kl  €  £.J( T /, /. ).  Note  that  the 
functions  A  kl  and  Ai;,  are  generally  different,  despite  the  fact  that  they  are  both  defined 
in  the  same  subset  T  k  ft  T  j.  of  R" 

It  is  clear  that  there  is  a  canonical  identification 


obtained  by  associating  with  A  =  (A )  g  LJ( Tj )  the  collection  (A/, )  defined  by 


^  K  |rh  L  =  •'/>  /.•  VI  €  Ik  . 
2 
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and  vice-versa. 

In  future  considerations,  it  will  be  convenient  to  use  both  notations  {X^l)  and  i^h)  for 
elements  of  £2(Fr).  No  confusion  should  arise  from  this  since  the  number  of  indices  used 
immediately  identifies  which  definition  of  L2{T j)  is  being  used. 

The  norm  of  I2(Fi)  is  the  natural  one.  namely 

/  \1/2  /  N  \t/2 

l^lo,rx  =  (  X]  )  ~  (  y.  I**' lu,r\rK  ) 

\k.L)€7  '  V=l  7 

We  shall  also  consider  the  spaces 

vK  =  {«-/,  € =o  in  rnr,.  |. 

with  norm  induced  by  the  norm  of  and  their  product 

Hl(P)  =  hnvh< 

equipped  with  the  product  norm,  denoted  by  ||  •  Ih  .n 

Elements  t'  of  H^(V)  can  be  viewed  as  AT  — tuples  (t \  )  with  17,  €  I  /,•  or  as  functions 
in  Q  (defined  almost  everywhere)  such  that  r/v  =  r|nfc.  €  V/v .  Both  points  of  view  will  be 
used  later,  and  once  again  no  confusion  should  arise  from  this. 

A  restriction  map  (trace)  is  defined: 


—  (vl<h.KU  )€  L2(Tj). 


which  is  evidently  continuous  for  the  norm  of  the  spaces  involved. 
There  is  a  canonical  embedding 


defined  by  v  1 — ►  (v^  )  where  .  On  the  other  hand,  the  space  £2(Ti)  can  be  split 

into  the  direct  sum 

iJ(rii  =  a+*a-. 


a+  =  idVi):r,,t  =rtA.  v(A'.i) g i) 


A_  =  {<rKi>:  rM.  =  -r,.K,  vfA'.ocr}. 

The  following  resalt  characterizes  the  element  of  Hj(fl)  among  the  elements  of 
Proposition  1.1:  Let  v  =  (t'j, )  €  Hq{V).  Then,  u  €  ffo(ft)  if  and  only  if  <  *'k,  X(. )  6  A+. 

Proof:  If  v  €  V(tt)  and  rr  =  n„  ,  it  is  obvious  that  v^.  =  17,  ,  i.e.  ((•/,.  )  € 

'  *  'rKL  ,r k  L  'ffc  Vf 

A+.  By  denseness  of  f2 )  into  H^iU).  continuity  of  the  trace  H?t(V)  Li{Ti)  and 
closeduess  of  A+  in  L2(Ti),  this  relation  extends  to  the  whole  space  #,',(0). 

Conversely,  we  have  Hq{V)  <—♦  L2(U)  in  t he*  obvious  way.  and  for  ^  6  P(ft)  and  t’  € 
ffi(P),  we  have 

v/here  =:  * '  * t'Kn)  i~  the  outward  unit  normal  vector  along  T Since  ^  s  0  on 

T,  we  find 


Y  f  =  Y  ( 

K=iJr«  a<.LmJr'iL 


and  the  terms  JpKl  cancel  out  front  the  hypothesis  vn  =  vl  on  IV/,  since  ///Vt  = 

- t'L , .  It  follows  that 

1. e.  dvfdxt  is  represented  by  the  Z2(fl)  function  defined  by  dv^  jdx,  in  1  <  A*  <  X. 

Thus,  v  €  HH ft),  and  since  r  =  0  in  0  T)  =  l\  we  have  v  €  //J  ( 0  >.  □ 

2.  The  space  A_. 

s 

Every  element  A  €  A_  can  be  viewed  as  an  element  of  ( HUT))'  =  II  \"!  via 

/x  =  I 


Vu’' = i  l 

h= i  ■,| » 


note  that 


(A, ti»)  =  0,  Vtr€#rJ(n).  VA  €  A_. 


T  /  ^KU'h  =  zl  I  ^}<U'hy 
K  =  1-/rfc\r  (K,D€X^rKt 


and  the  terms  /P|fi  A^u'/c  and  frKL  A/.u'i,  cancel  out  because  ir/v-  =  ti'i,  on  T/vl  (see 
Proposition  1.1)  and  A^-  =  —A/,  on  Tkl  if  A  €  A_. 


Definition  2.1:  The  space  A-  is  the  closure  of  A_  in  =  II  Vfc. 


By  definition,  every  A  €  A_  is  the  limit  in  {Hq{V)Y  of  a  sequence  A;  €  A-.  Since 
(A*,tr)  =  0,  Vu*  €  Hq(Q)  from  (2.1).  we  obtain  by  continuity  that 


(A,  «*)  =  0,  VA  €  A_.  Vu><=  #<!(*?). 


Let  i4  —  (a,j)  6  (Lao(Q))n><n  and  let  Hq  A(Q)  be  the  subspace 


<„(!!)  =  !>•  €  H‘m :  V  ■  .-tVr  e  I2(n)}, 


equipped  with  the  norm 


MI.M.n=(IMI?.n  +  |V--4V1.|5.„)|/J. 


Proposition  2.1:  The  mapping 


.•€  P<ni. 


m- 


(-4Vi  /,  i',, )  6  A- 


(where  t =  i'|„  ,  1  <  A"  <  .V)  can  be  extended  as  a  linear  continuous  mapping  from 
into  A_ . 


•  i 


A 


Proof:  It  is  plain  that  [di’K /dvA )  e  A-  as  soon  as  i>  €  P(  Sf ) .  In  addition,  by  Green's 

formula. 

/  T^'fk  =  /  V.4,-k  VwK  +  /  V  (.4Vra  V.rK  €  l'A 
ir*  ^  JUh  Jnh 

Thus,  adding  up  these  relations,  we  find 

1  ((fe)’U,)  n  Mi.*.  Vu’  €  H^V)- 

where  C  >  0  is  a  constant  depending  only  upon  .4.  Hence, 

ll(|^-)lkr  <C|M|,-4,o,  VretHO). 

where  ||  •  ]|.  r  denote  the  norm  in  {H^(V)Y.  That  (dv^/dt'A)  ran  he  extended  as  an 
element  of  A-  for  v  6  A{U)  follows  from  this  inequality  and  denseness  of  ViCl)  into 
//J ^(fi)  (for  the  ||  •  HiM.fl  norm)  by  standard  arguments.  □ 

,v 

Remark  2.1:  Because  A-  C  (Hq{P)Y  =  II  IT.  elements  of  A.  are  *V- tuples  A  =  (A/v ) 

l\  =  i 

with  A^  €  VL-  but  unlike  elements  of  A_,  they  cannot  be  identified  with  collections  (A  /%-  /, ) 
for  (I\,L)  €  I  (since  elements  of  VA  are  not  functions,  they  cannot  he  restricted  to 
subsets).  □ 


3.  Application  to  error  estimation  in  linear  problems. 

In  addition  to  the  coefficients  .4  =  of  the  previous  section,  we  also 

consider  6  =  (6.)  €  and  r  €  and  set 


a(r.  u’ ) 


=  /  AVr  ■  Vie  4-  f  (4  Vv)w  +  /  cvw,  Vr.  w  6 

J  fl  J  n  J( I 


We  shall  henceforth  assume  that  «(■.■)  above  satisfies  the  conditions  (1  la/h)  Among 
other  things,  this  ensures  that  given  /  e  L2(tt).  there  is  a  unique  solution  u  €  ( fi i  of 

the  problem 


n(u.  >') 


Vue  Him. 


Observe  that  u  solves  the  PDE 


(3.1) 


— V  ■  (rlVu)  +  fc  Vu  +  cu  =  f  in  0. 


whence  V  (eiVu)  =  b  Vt /  +  cu  -  f  €  £J(fl ).  i.e.  ti  €  //,j  ^(0). 

For  1  <  A’  <  N,  let  '  1  k  x  l  a  — *  R  be  the  (contiimons)  bilinear  form 

“A(t'A.U’A)=  /  AVl't,  Vll'A  +  I  (A  Vt'Altl’A  +  [  CI'AII  A-  Vt'A.WA  €  l  A 

lit*  In,  Vnk 

We  assume  that  each  bilinear  form  a  a  satisfies  the  conditions  (l.la/b),  so  that  given 
u  €  tf,j(fi)  the  problem 

(3  2)  Otdfe.t'ti  )  =  /  /t’K  =  OA(ti.t')  +  (Ak.I'a).  Vo,,  €  I’a 

•/  n, 

has  a  unique  solution  4a(Aa)  €  1  a  for  every  A  a  6  V,'.  In  particular,  given  A  =  (Aa)  6 
A-,  there  is  a  unique  0(  A I  €  //0'(P)  such  that  0yv  =  0  y,(  A  a)  solves  (3.2)  for  1  <  A’  <  .V. 
The  following  theorem  shows  that  the  error  t/  —  ti  is  characterized  by  a  property  of  0(A) 
holding  for  one  and  only  one  A  6  A  _ . 


Theorem  3.1:  The  following  conditions  are  equivalent : 

(i) *(A)€  ffj(«). 

(ii)  ( <i>  a(  A  a  )|rK  r )  €  A+, 

(iii)  0(A)  =  u  -  ti. 

(iv)  A  =  (duK/dvA). 

Proof:  For  (i)  •»  (ii).  see  Proposition  1.1.  Before  proving  that  (i)  <=>  (iii)  (iv).  let  us 
note  that  from  (3.1)  and  Green's  formula,  we  have 

/  /"A  =  Oa('M'a)-  ( ’  t'A  )  • 

jo* 

where  du^/dvA  G  l ^  vanishes  in  f/J ( Q y, )  (the  validity  of  Green’s  formula  above  is  due 
to  u  €  Hq  By  substitution  into  (3  2).  we  get 

(3  3)  <ia(^a(  Aa  )  t'y,  )  =  oy,(ti  -  ti,  i'a)  +  (Aa  -  -5— fy, ).  Via  €  l  a. 

w, 


and  hence 


N  N  Ndu 

)  =  ^ok(u  -  u. rjv)  +  ^(A/,  -  y^-.t'h).  Vr  e  H,!(P) 

h=l  h*|  h=|  '* 

(i)  ^  (iii):  Choose  t»  €  Hq  ( (2 )  in  the  above  relation.  From  (2.2)  and  A.  (du^  fdv\ )  €  A_, 
it  follows  that 

N  N 

=  y  nh ( 11  -  «•>•/>)•  Ve  e 

A=l  />  =  ! 

But  u  -  ti  6  //q ( 0 )  and  0(A)  6  Hi(Q)  from  (i).  so  that  the  above  equality  roads 
a(0{A).r)  =  a(u  -  u.f).  Ve  €  ff, J(fl). 
i.e.  0(A)  =  u  -  u  since  a(  .  )  satisfies  ( l.la/b). 

(iii)  =>  (iv):  If  d>( A )  =  u  -  u,  tlien  0h(A/,  )  =  (1/  -  i))|n  and  (3.3)  yields 

(As  -  ~.  ra )  =  0.  Via  els.  1  <  /t  <  .V. 

O'-.j 

As  a  result.  Ah  =  du^/'di >a  (equality  in  ),  1  <  l\  <  X ,  i.e  A  =  (du^/d i'a). 

(iv)  =>  (i).  If  A  =  (du^/d i/aI,  then  Ah  =  3uh/9«'*,l  <  K  <  .V,  and  (3.3)  shows  that 
*A'Oli )  is  charactcrizeed  l>y 

Oh(^a(Ah  ).i'h  )  =  «h(u  -  ii.e/v  ),  Vi'h  €  V), 

But  certainly  0;v(  A/,  )=■(«-  m  )|n  satisfies  the  above  relation,  and  lienee  is  the  unique 
solution  of  (3.3)  This  shows  that  0(A)  =  11  -  u  c  W,J($1)  CJ 

Corollary  3.1:  The  functional 

A  €  A-  • — •  7(A)  =  -  y  kh(A;,  )  --  */.(  A/,  )|Jj-Kl 


has  the  unique  minimizer  A  =  {du^- JOi'a).  for  which  <2»( A )  —  u  -  u. 


Proof:  From  (iv)  ^  (ii)  in  Theorem  3.1.  we  have  J{{du  ^ /Oi'a))  ~  0.  and  hence  A  = 
(dua/dt'*)  is  a  lanumizer  of  7  since  7  >  0.  This  also  shows  that  7(A)  =  0  for  any 
minimizer  of  7  but  then  A  —  (du^/di'A)  from  (ii)  =$>  (iv)  in  Theorem  3.1.  Thus.  A  = 
{duh  is  the  unique  minimizer  of  7.  and  0(A)  =  t<  -  u  by  (iv)  =>  (iii)  in  Theorem 

3.1  □ 

Corollary  3  1  suggests  a  strategy  to  calculate  the  error  ii  ».  hy  minimizing  tin*  func 
tional  7  in  (3.4).  It  is  important  to  notice  that  the  functional  J  is  quadratic.  More 
specifically,  denote  by  7*  €  \\  the  (unique)  solution  of 

.»**)  =  /  /tK  -  aglN't'A  ).  Vr/cClV- 

7uk 

and  let  €  C{V^  ,\'h  )  be  defined  hy 

a#\ ( I  k rK •  *’#»  1  =  <f,/vr/0-  v‘‘/x  €  '/!•  Vi/x  6  Vh 
Then.  Ok  ( Ag  )  —  I’a  A/x  -f  7/x .  1  <  A*  <  A’,  for  A  c  A-.  and 

^  5Z  { K  h Ak  -  t  | , A|^ Iq  jKt  +  2  f  (t\A;x  -  I  /  A/ )(7k  -  7;  ) -f  |7k  -  7/,|J  , 

"tK.tjtr  1  *'r*t 

k  >/. 


that  is. 


13.51  Zl\|-  i»(VAK  /  (['k*!, ,1.^  -  .1,  )+i  Y.  M».  -  M'r.,- 

h>l. 


where,  for  A./j  C  A-  we  have  set 


*(  A. /i)  =  E  /  /x  Mv  -  1 7  A/  mI’k /;/v  I  i.  M  ) 

(  h.Dfl  *  L 
h>i 

n 


I 


Lemma  3.1:  The  bilinear  mapping  a  is  an  inner  product  in  A-. 

Proof:  It  is  obvious  that  o  is  bilinear  and  symmetric.  That  <r{  A,  A)  =  0  =>  A  =  0  follows 
from  (ii)  =>  (in)  in  Theorem  3.1  for  the  special  case  when  /  =  0  (hence  u  =  0)  and  u  =  0 
in  (3.2).  Q 

The  norm  in  A_  induced  by  the  inner  product  <»(,•)  will  henceforth  he  referred  to  as 
the  “<r-norm”.  It  is  easily  seen  that  A-  cannot  be  complete  for  the  <r-norm.  and  lienee  wfe 
shall  call  A*  the  completion  of  A  for  the  (* norm.  From  the  trace  theorems,  there  is  a 
constant  C  >  0  such  that 


0(.V.M,/2  <  C  MUr.  VX€ 
i.e.  the  it -norm  is  weaker  than  the  ||  •  ||.  p  norm  in  A-. 

Theorem  3.2:  The  functional  J  is  continuous  and  coercive  in  A_  equipped  with  the 
(r*iuirtn.  As  a  result,  it  can  be  extended  as  a  continuous  and  coercive  quadratic  functional 
in  A®  with  the  same  minimum  value  (that  is.  0).  In  particular.  A  =  [duii/di'.\)  remains 
the  unique  minimizer  of  .7  in  A* 

Proof:  Continuity  and  roereivity  c»f  .7  in  A-  for  the  a  norm  is  trivial  from  (3.5).  That 
it  can  be  extended  to  At  with  the  same  minimum  value  is  also  trivial  (using  denseness  of 
A-  in  its  completion  At)  It  thus  remains  uounegative  and  has  a  unique  minimizer  in  At, 
and  this  minimizer  must  be  since  J({du/S  /du^))  —  0.  □ 

We  now  face  the  key  problem  in  following  the  alxive  approach  In  practice,  the  minimizer 
[duiy  /Oi'A )  of  J  will  be  obtained  via  minimizing  sequences  A;  6  At.  By  denseness  of  A_ 
in  At.  such  sequences  can  be  chosen  in  A-  (and  even  in  A-  by  denseness  of  A_  in  A_). 
ami  for  each  index  j  we  obtain  an  element  6{\})  £  Hq[P).  It  is  easily  seen  that  0{X}) 
need  not  approach  u  -  i«  in  H^i  P)  as  j  tends  to  x.  This  is  related  to  the  fact  that  the 
<i  norm  is  t f m >  weak  to  ensure  continuity  of  the  bilinear  form 

N 

ia.i  ) e  a.  -  n’tr i  —  £(\h.rh) 

i- =, 

10 


» 


» 


» 


» 


m 


when  A_  is  equipped  with  the  e-norm  Thus,  in  practice,  minimizing  tlic  functional  J  using 
minimizing  sequences  (the  only  possible  option  snce  (dux /dvA)  is  of  course  unknown)  will 
not  necessarily  provide  the  desired  result,  i.e.  will  not  yield  an  approximation  on  u  —  u  in 
Hq(V).  Fortunately,  this  difficulty  can  be  overcome  because  of  the  following  result. 

Theorem  3.3:  Suppose  that  u  £  for  some  s  >  3/2.  Then. 

(i)  There  are  minimizing  sequences  (A;)  of  J  such  that  \J  £  A_  and  ( AJ)  is  bounded  in 
A.,  (that  is,  bounded  in  £J(Fj)). 

(ii)  If  (AJ)  is  a  minimizing  sequence  of  J  such  that  AJ  £  A-  and  { A^ )  is  bounded  in  A-, 
we  have 

lim  ||0(A;)  -  (»i  -  u ) ||  i  r  =  0 

)—oc 

Proof:  If  u  €  H‘( fi)  for  s  >  3/2.  then  dmJdvA  e  H’~\(  )  C  A2(Ta  ),  1  <  I<  < 

N.  Thus,  in  particular,  (3u/,  ldvA)  6  L2(rr).  Also,  if  u  6  C^ffi),  it  is  obvious  that 

(3uk /dvA )  €  A_,  and  this  relation  remains  valid  for  u  €  H‘(Q)(s  >  3/2)  by  denseness  of 

C*(ft)  into  H‘( fl).  continuity  of  the  mapping  v  6  H‘(Q)  • — >  (dv^/d vA)  €  I2(T j)  and 

elosedness  of  A.  in  L2( Tj )  Existence  of  minimizing  sequences  (V)  with  At  €  A-  and 

(AJ)  bounded  in  A-  follows  at  once  from  this  result  (choose  e.g.  At  —  (du^ jdvA ).  V)). 

Next,  let  ( AJ )  be  any  minimizing  sequence  of  J  such  that  At  €  A  _  and  ( AJ )  is  bounded  in 

N 

A-.  By  compactness  of  the  embedding  I2(r;v  \  T)  IV  ,  there  is  A  €  (/fd(P))' =  II  17, 

K*1 

and  a  subsequence  (A'*)  such  that  lim  || AJ*  -  A||.  r  =  0.  Thus.  A  €  A-  by  definition 

A  — oo 

of  the  space  A*.  Since  the  functional  J  is  continuous  for  the  <7-norm  (see  Theorem  3.2). 

hence  for  the  stronger  ||  ||.  p  norm,  we  have  J{  A)  =  lim  J{  At* ),  so  that  A  minimizes  J. 

A  — -X, 

It  follows  that  A  =  ( du /,  /dvA )  (see  Theorem  3.2). 

The  above  rev,  als  that  ( du /drA )  is  the  unique  cluster  point  of  the  sequence  (A;)  in 
A-.  SO  that  the  whole  sequence  (At)  converges  to  (dti^/di'A)  in  A..  This  amounts  to 
say;ng  that  \’K  tends  to  du^/d i'A  in  \"h  for  every  1  <  A"  <  .V  Now.  from  (3.3)  we  have 

nh  ( )  -  ('<  -  u).  f#,  )  <  HA},  -  ^j“ll*.nK  ll'  l.  III  Ok  •  v,'/>  e  1 /,  ■ 


Taking  the  supremum  over  vk  €  \’k ,  II<’K  II i  Hie  <  1 .  and  using  the  fact  that  (!/,(■,)  satisfies 
the  conditions  (l.la/b),  we  get 

II*k(A{J  -  (u  -  »)||i.(tK  <  Ch||A'K  - 

where  Ck  >  0  is  a  constant  independent  of  j.  Thus, 

(TMM  -  (U  -  »)li;,ne)  - 

with  C  —  max  C*  ,  and  the  right-hand  side  is  just  C||AJ  —  [  J  ||,j»  and  hence  tends 

\<h'<N  \  A  / 

to  0  as  j  tends  to  oc  as  noted  earlier  in  the  proof.  D 

By  a  method  of  proof  similar  to  that  of  Thoerem  3.3.  we  obtain  the  following  result 
without  any  extra  assumption  of  regularity  for  the  solution  u: 

Theorem  3.4:  0)  There  are  minimizing  sequences  (A;)  of  J  with  AJ  €  A-  and  (A>) 
bounded  in  A_ . 

(ii)  If  (A>)  is  a  minimizing  sequence  of  J  such  that  \}  6  ft  and  (AJ)  is  bounded  in  ft, 
we  have 

lim  ||^(  AJ)  -  (u  -  u)||r.r>  =  0. 

)—oo 

.V 

for  every  0  <  r  <  1,  where  ||  •  ||r,r  is  the  product  norm  of  the  space  Hr(V)  =  ^11  Hr(Qk  ). 

Proof:  (i)  is  trivial  by  considering  the  constant  sequence  (dux ('»)  (A->)  be  a 

minimizing  sequence  in  A-  which  is  bounded  in  A_,  and  let  A  be  a  cluster  point  of  (A;)  for 

N 

the  weak  topology  of  A_,  ( A_  is  a  closed  subspace  of  the  Hilbert  space  If  (A;k) 

i**  h  subsequence  tending  to  A,  we  still  have  J(A)  =  lim  J{  AJfc )  =  0  because  J  is  convex 

t— oo 

and  continuous  for  the  ||*||.p  :’orm,  hence  weakly  lower  semicontinuous.  Thus,  once  again. 
X  -  {dux/du/i)  and  the  whole  sequence  (A;)  tends  weakly  to  (du  ^  /duA).  At  this  point, 
note  that  together  with  this  result,  (3  3)  implies  that  0k  (A^  )  tends  weakly  to  (u  —  u)|OK 
m  Ik  If  £ h  €  £»£(!*)  is  defined  by 

afcft’K.U'K  )  =  (LhWK'Wk ),  VPK-,tt'A‘  €  I  A  . 
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where  (■,■)  denotes  the  inner  product  of  H'(Qh).  we  have  <t>h(  )  —  (u  —  u)|„k,  =  L  a 

where  p’h  €  Vj,  is  defined  by 

(?/,,'/> )  =  (K  ~  ••/.),  v,'a  £  '  a 

That  ffJh  — *  0  in  IV  follows  from  A£  —  du^ /dv,\  — *■  0  in  \r,h  and  hence  <£h(AJ )  — ■  (u  —  w)|nfc, 
by  continuity  of  L^  . 

The  theorem  now  follows  from  compactness  of  the  embedding  Hr(Uh)  *— *  Hl(Qx)  for 

0  <r  <  1,  1  <  K  <  N.  □ 

The  practical  aspects  of  Theorems  3.3  and  3.4  is  captured  by  the  following  corollary. 

Corollary  3.2:  Let  B  C  A-  be  any  closed  ball  such  that  (du^/d v*)  €  B.  Then.  B  is 

closed  in  Ka  ,  inf  J{A)  =  inf  J{  A)  and  if  \*  £  B  is  a  minimizing  sequence  of  J\B ,  we  have 
A€B  A€A* 

lim  ||<£(A;)  -  (u  -  n )|| r.p  =  0, 

J—OO 

for  every  0  <  r  <  1. 

Furthermore,  if  u  €  for  some  s  >  3/2  and  B  C  A_  is  any  dosed  ball  such 

that  (duh/duA)  €  B ,  then  B  is  dosed  in  Al,  infJ(A)  =  inf  J(A)  and  if  A;  6  B  is  a 

AeB  A€A1 

minimizing  sequence  of  Jifl .  we  have 


lim  ||tf(  AJ)  -  (it  -  t'OHi.p  =  0. 


Proof:  That  B  is  closed  iu  Al  follows  from  weak  compactness  of  B  in  its  ambieut  space: 
If  A  *  €  B  and  \J  — *  A  £  A"  in  A*,  let  //  €  B  be  such  that  \)k  — *  //  in  A-  or  A-.  In  both 
cases,  this  implies  \lk  — *  /t  in  A"  sine-'  the  embeddings  A_  ‘—♦A  <—*  Al  are  continuous, 

whence  /i  =  A  and  A  6  B.  That  the  infimum  of  J  is  the  same  in  B  and  in  the  whole  space 
Al  is  trivial  from  {du^/du)  6  B  being  the  unique  minimizer  of  J  in  Al.  Convergence  of 
<p{\})  to  u  —  u  follows  from  Theorems  3.3  and  3.4.  □ 

In  practical  applications,  the  hypothesis  that  some  ball  B  containing  (duix  /Bi'a)  is 
known  is  not  a  severe  limitation  since  some  useful  information  can  be  derived  from  the 
approximate  solution  ».  Sec  Section  4. 


13 


w 


t a 


4.  Practical  Aspects. 

In  practice,  the  appoximate  solution  u  will  be  a  finite  element  approximation  ti*  6 
V'*  where  V  *  is  a  finite  element  subspace  of  Hq ( Q )  and  the  partition  V  corresponds  to 
the  partition  of  0  into  elements.  It  is  certainly  not  restrictive  to  assume  that  u*  is  a 
“reasonable  approximation  of  u.  as  opposed  to  a  completely  random  element  of  V'*,  and 
hence  that  despite  the  fact  that  ||t<  -  u* ||,.0  may  not  be  small,  or  small  enough,  still 
the  derivative  (dukk/dvA)  gives  some  idea  of  where  in  A-  (or  A-  if  u  6  a  > 

3/2)(chig-  ldvA )  is  One  problem  is  that,  in  general,  {dit)jdvA )  $  A_,  but  this  is  easily- 
remedied  by  replacing  (du‘  /dt'A)  by  say,  (§“-)1/2  =  ({<dukk/dvA  -  dhJdVA)\K, t)eX. 
Other  weighted  averages,  as  suggested  in  [  }  ran  also  be  used. 

The  ball  B  of  Corollary  3.2  can  be  chosen  as  any  ball  with  center  (du/c h'a)i/2  and 
arbitrary  radius  B  >  0.  Because  of  denseness  of  A-  in  A_  and  A*  ,  finite  dimensional  ap¬ 
proximations  of  the  space  At  can  be  obtained  by  choosing  a  scale  Ai  of  finite-dimensional 
subspaces  of  A-.  One  possible  choice  is  given  by  collections  A  =  (A;,)  such  that  A;,  is  a 
polynomial  of  degree  <  j  on  each  face  of  V, k  not  lying  on  the  boundary  of  fl.  It  should  be 
kept  in  mind  that  membership  to  A_  requires  the  condition  A =  -A/  .  and  hence 

rKL  1 rKL 

A  a  and  XL  cannot  be  chosen  independently.  Because  no  continuity  condition  is  required  of 
A;,  at  the  intersection  of  two  faces,  bases  for  such  spaces  are  easily  found.  For  instance,  if 
}  =  0.  a  basis  is  given  by  the  collections  (Am)  where  XKL  =  0  if  (A'.Z)  /  (A'0.l0)  where 
(ho.Lo)  €  I  satisfies  A  a  >  Lo,  A  =  1-  -  -1.  letting  ( K0.  L0  )  run  over  all  such 

pairs.  It  should  be  clear  how  this  procedure  can  be  extended  to  obtain  bases  for  arbitrary 
polynomial  degree  j . 

One  may  then  define  X1  to  be  the  minimizer  of  J  in  HO  Ai,  a  closed  convex  subset  of 
Ai  The  sequence  (At)  is  a  minimizing  sequence  for  J  in  B.  and  hence  the  conclusion  of 
Corollary  3.2  regarding  convergence  of  <t>(  X> )  to  u  -  u*(  =  u  -  ti  here)  applies  Rather  than 
minimizing  J  in  D  n  Ai  -  0  a  constrained  problem  -  it  seems  advisable  to  minimize  J 
in  .V  (a  linear  problem)  and  check  whether  this  minimizer  is  in  B.  and  modify  the  choice 
of  the  minimizer  only  if  this  is  not  the  case  (i.e.  the  minimizer  is  unreasonably  far  from 
{duk/d 

Naturally,  replacing  A-  by  a  finite  dimensional  approximation  .Vi  is  not  enough,  since 


H 


0(A)  cannot  be  calculated  exactly.  However,  since  0(A)  =  (0*  { A*  ))  is  obtained  by  solving 
independent  local  problems,  numerical  approximations  can  be  obtained  at  a  low  cost  with 
high  accuracy.  In  other  words,  the  error  between  0(A)  and  its  numerical  approximation  may 
be  considered  negligible  when  compared  with  ||u  —  «*||i and  hence  the  convergence  result 
of  Corollary  3.2  may  be  used  safely  with  0(A)  replaced  by  its  numerical  approximation. 

Remark  4.1:  When  the  spaces  Ai  are  those  previously  described,  i.e.  their  elements  are 
collections  of  polynomials  with  degree  <  j  on  the  interfaces,  it  is  desirable  that  the  space 
chosen  for  the  approximation  of  0*(A *  )  when  A  =  (A*-)  €  Ai,  contains  the  polynomials 
of  degree  <  j  +  1  in  ft*.  Indeed,  A*  is  used  to  obtain  an  approximation  of  du^/d 
whereas  0  *(  A  * )  represents  an  approximation  of  (u  -  iifc)(nfc,  •  an(*  hence  polynomials  with 
degree  at  least  j  +  1  should  be  used  to  approximate  0*(A*)  for  consistency  (assuming 
A  =  (ai;)  essentially  constant  ou  each  element).  The  above  considerations  also  show  that 
j  +  1  should  be  at  least  equal  to  the  degree  of  the  polynomials  used  for  the  calculation  of 
which  comes  as  no  surprise.  □ 

5.  Nonlinear  problems. 

Let  a{  ,  ),a*  (  ,  )  be  the  same  bilinear  forms  as  in  Section  3,  and  let  F  :  ft  x  R  — ►  R 
be  a  mapping  such  that  F(\  )  is  of  class  C1  for  almost  all  x  €  ft  and  F{-,y)  is  measurable 
for  every  y  €  R.  We  also  assume  that  there  are  constants  {3  >  0,  >0  such  that 

|0,C(T.y)|  <  d  +  7|y|p'2 

for  almost  all  x  €  11  and  •■very  y  €  R.  where  p  >  2  is  a  real  number  satisfying  p  <  jf 
n  >  2.  This  implies  that  the  mapping 

F  :  u  €  Lp{,U)  =>  F(u)  €  L”  (fl)  (p*  =  p/(p  -  1|) 

defined  by  F(  w )( -r )  =  F(j,ti(x))  is  of  class  C  with  derivative 

DF{u)h  =  DtF(x.u(x))/i(x). 

IS 


ft 


ft 


ft 


ft 


ft 


ft 


ft 


It  then  follows  from  Taylor’s  formula  that 


|F(u)  -  F( t')|o.,..n  <  C((mcas  Q)1?  +  7(l<,2n  +  |t-IS*n)l«  -  «’lo.,.n. 


for  every  pair  (u,i»)  €  (LP(Q))2,  where  C  >  0  is  a  constant  independent  of  u  and  v. 

This  shows  that  if  uh  €  V*  is  a  finite  element  approximation  of  u  6  such  that 


lim||u  -  tr* ||i ,n  =  0, 


..  |u  ~  u*lo.P,n 


=  0, 


we  have 


(5.1) 


|F(u)-f’(u‘)lo.,.,o  <€{h)  || u  -  u‘||i.n. 


with  lime(h)  —  0 

h  —0 

Suppose  now  that  u  solves  the  problem 


o(u,t’)  + 


f  F( ..).>  =  / 
J  a  J  n 


/<’ 


i.e.  u  €  Hq{Q)  and 


-V  •  (/IVu)  +  6  ■  Vu  +  cit  +  F(u)  =  f  in  fl. 


In  analogy  with  what  was  done  in  Section  3,  let  A  €  A-  and  solve,  for  1  <  K  <  iV. 


I  aft  (dafAa  l.t'/,  ) 
1  VeK  e  V 


-  «!#>•(  11* 


'•/> ) 


/la 


F(tih)rh  +  (Aa.t’/,) 


that  is. 

(  aa  (0a  (A/,  ),t’a  )  =  a;v(ti  -  ti*,t'a  )  +  /  (F(u)  -  f’(u*))i7>  +  (A/,  -  -z^-.vk) 
(5.2)  f  JaK  Ova 

\  Vi  a  €  V 
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By  the  method  of  proof  of  Theorem  3.1,  it  is  easily  seen  that  the  only  case  when  0{A)  = 
€  Hq(Q)  is  when  <£(A)  coincides  with  the  solution  t/'*  €  /fg(Q)  (not  Vh)  of 


> 


(53)  a(V*,v)  =  o((u-u‘ ),«•)+  f  (F(u)~  F(uk))v.  Vt  >e//<j(n), 

Ju 

and  that  if  so,  A/,  =  3(t/’*  +  uK)^/dvA.  1  <  K  <  N.  In  fact,  in  terms  of  tc*  above.  (5.2) 
reads 

I  0<v(^A  (^a).  I’/,)  =  «a(«'\  t’A  )  f  (Af, - - ~'vf) 

\  VfA  6  IV 

This  relation  can  be  used  to  show  that  if  (AJ)  is  a  minimizing  sequence  of  the  functional  J 
of  Section  3  (which  now  has  unique  minimizer  (d(4’h  +  u>‘)ix/duA ))  and  if  (A>)  is  bounded 
in  A_,  then 

(5.4)  lim  ||^(A'’)  -  0*||i,p  =  0. 

J— oo 

(And  if  (A;)  is  only  bounded  in  A_, 

lim  ||d(AJ)  -  A'*||rT-  =  0  ,  V0<r<l.) 

)—  00 

To  obtain  an  estimate  for  ||0(AJ)  -  («  -  u*}||,  Vy  write 

(5.5)  IWAJ)  -  (u  -  u*)||,,,.  <  ||*(A')  -  v*||i.p  +  |l( tz  -  u*)  -  v*l|i,n. 


> 


I 


From  (5  3), 

a(*k  -  (u  -  u*).r)  =  ( (Flu)  -  F(ul'))v. 

Ja 

Taking  the  supremum  over  r  €  ||c||iii  <  1  and  using  the  fact  that  <i(  .  ■)  satisfies  ® 

(l.la/b),  we  get 

<%>*  -  ( u  -  u*)||i,n  <  sup  f  (F(u)  -  F(<il'))r, 

II ell i  n<l  Ja 


where  c  >  0  is  a  constant  independent  of  u  and  h.  Next,  using  |jjl(/'(u)  --  F(u*))t>|  < 
(F(u)  -  F(uA  )|0y  .nWlo.p.n  ®^d  continuity  of  the  embedding  <—»  Lp(ft),  we  obtain 

||</.‘-(«-u‘)||..Ii  <C|F(«)-F(U‘)|.r,fl, 
where  C  >  0  is  a  constant  independent  of  ti  and  h.  Finally,  using  (5.1)  we  see  that 

111/'*  -  (u  -  u*)||,,n  <  c<(h)||u  -  u‘||i.n- 
Substituting  into  (5.5),  we  get 

IMA')  -  (u  -  u*)||,,*  <  ||i/'(A')  -  V'*||.,r  +  CV(A)||U  -  U*||,.„. 
and.  in  particular,  if  (5.4)  holds: 

(5.6)  lim  IMA')  -  (u  -  u*)||,,j>  <  Cr(ft)||u  -  ti*||,.n, 

)~>oo 

i.e.  || c»( )  —  (u  —  u*)||]i7  becomes  negligible  with  respect  to  ||u  —  up||j  q  when  0(A->)  is  a 

good  enough  approximation  of  0*  because  lime(/i)  =  0. 

h  — o 

Remark  5.1:  From  (5.6),  it  also  follows  that  uh  +  0(AJ)  is  a  better  (enhanced)  approx¬ 
imation  of  u  than  uh .  Interestingly,  this  enhanced  solution  is  obtained  by  solving  only 
linear  problems  (calculation  of  0(A->)).  □ 
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